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Spin pumping from a ferromagnet into a hopping insulator: the role of resonant 

absorption of magnons 
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Motivated by recent experiments on spin pumping from a ferromagnet into organic materials in 
which the charge transport is due to hopping, we study theoretically the generation and propagation 
of spin current in a hopping insulator. Unlike metals, the spin polarization at the boundary with 
ferromagnet is created as a result of magnon absorption within pairs of localized states and it spreads 
following the current-currying resistor network (although the charge current is absent). We consider 
a classic resonant mechanism of the ac absorption in insulators and adapt it to the absorption of 
magnons. A strong enhancement of pumping efficiency is predicted when the Zeeman splitting of 
the localized states in external magnetic field is equal to the frequency of ferromagnetic resonance. 

Under this condition the absorption of a magnon takes place within individual sites. 

PACS numbers: 85.75.-d,72.25.Rb, 78.47.-p 


I. INTRODUCTION 


The phenomenon of spin pumping from a ferromagnet 
(F) into a normal (N) layer is one of the most promi¬ 
nent approaches to the generation of pure spin currents. 
A prime manifestation that pumping indeed takes place 
in realistic F-N structures is the additional broadening 1 
of the ferromagnetic resonance (FMR) in F, caused by 
a contact with N-layer. This additional broadening 
was first observed experimentally in Ref. [2]. Another, 
more delicate, manifestation of pumping was reported 
shortly after. Namely, the injected spin current, entering 
the nonmagnetic material with spin-orbit coupling (like 
Pt) causes a voltage drop across the current direction. 
This voltage drop is due to the inverse spin-Hall effect 3 
(ISHE), and has a maximum when the frequency of the 
microwave radiation driving the ferromagnet, c o, is equal 
to the FMR frequency, uj fmr . Pioneering observations of 
pumping via ISHE in Refs. [4-6] utilized Pt as the nor¬ 
mal layer. 7 9 They were followed by reports on similar 
observations of pumping into different materials 10-13 , in¬ 
cluding prominent semiconductors GaAs 14 , Si 15,16 , Ge 17 , 
and, most recently, graphene. 18 Experimental results on 
the electric field generated due to ISHE , E ishe , are an¬ 
alyzed using the relation E ishe oc J^ x cr, where J(S) 
determines the spatial direction of the spin current flow 
and its magnitude, while cr is its polarization. The mag¬ 
nitude of the spin current is given by 


J(S) 


g n C m(t) x 


dm(t)~ 
dt - P 


(i) 


where z axis is taken along the static part of the mag¬ 
netization. In Eq. (1) the constant C characterizes the 
properties of the normal layer (like ratio of thickness to 
the spin-diffusion length), while m(t) describes the mag¬ 
netization dynamics in the ferromagnet. The expression 
for J has the same form as the damping term in the 
equation that governs m(t). It was a remarkable experi¬ 
mental finding 5 that ISHE voltage exhibits essentially the 
same behavior as a function of microwave power and the 


deviation of uj from uj FM r as the additional FMR damp¬ 
ing. 

Microscopic physics of pumping is encoded in the mix¬ 
ing constant 1,19,20 in Eq. (1). A fundamental process 
underlying the pumping is the inelastic electron-magnon 
scattering at the F-N interface. Microscopic treatment 
of this scattering 21,22 assumes that electrons of the nor¬ 
mal layer impinging on the interface with ferromagnet are 
plane waves. On the other hand, in a number of recent 
papers 23 26 spin pumping into organic materials sand¬ 
wiched between ferromagnet and Pt has been reported. 
Strong temperature dependence of the resistance in these 
materials 27 suggests that the charge transport is due to 
hopping of polarons 24,26 , so that the description of pump¬ 
ing based on plane waves does not apply. This raises the 
question about the microscopies of spin pumping in the 
localized regime. 

In the present paper we consider theoretically the spin 
pumping into a hopping insulator using the minimal 
model of coupling of localized states to a ferromagnet. 
We demonstrate that, unlike metals, the underlying pro¬ 
cess responsible for pumping is the resonant magnon ab¬ 
sorption accompanied by transitions between localized 
states, see Fig. 1. A distinctive feature of pumping into 
an insulator is that that the pumping efficiency, com¬ 
monly described by a constant, depends strongly on 
the external dc magnetic field. This is because, in addi¬ 
tion to causing the spin precession in ferromagnet, this 
field modifies the spin structure of the localized states 
between which the magnon is absorbed, see Fig. 1. The 
effect of external field is most pronounced when the wait¬ 
ing time for a hop is longer than the period of the ac field 
which drives the FMR. Since the resonance frequency, 
uj fmr , depends on the orientation of the external field 28 , 
for certain orientations 29 this frequency coincides with 
the Zeeman splitting of the localized states, Fig. 2. Spin 
pumping is most efficient for such orientations, since the 
absorption of magnon takes place within individual sites. 
We also show that, with no charge current, the spin po¬ 
larization generated at the F-N boundary, spreads in the 
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FIG. 1. (Color online) Elementary processes underlying the 
spin pumping into a metal (a), and into an insulator (b). In 
the metal, an f electron, impinging on the N-F boundary, is 
primarily reflected elastically with amplitude rq. Spin pre¬ 
cession in F gives rise to inelastic reflection with amplitude 
associated with the emission of a magnon. A f electron 
is either reflected elastically with amplitude rq, or inelasti- 
cally, after absorbing of a magnon, with amplitude f|q. The 
injected spin current is proportional to In the in¬ 

sulator, only inelastic processes are at work. Emission and 
absorption of magnons take place within pairs of localized 
states. 


insulator along the same percolation network 30,31 that 
determines the electrical resistance. 


II. ABSORPTION OF MAGNONS AT F-N 
BOUNDARY 



FIG. 2. (Color online) Illustration of pumping in metal (a) 
and in insulator (b) in the presence of a Zeeman splitting, 
A z . In metal, the absorption (emission) of a magnon, huo , 
near the F-N boundary does not conserve momentum, and 
thus is insensitive to the ratio A z /hu. By contrast, in and 
insulator, and near the condition hw = A z , the absorption 
(emission) of a magnon is resonant 


magnon absorption takes place between the discrete lev¬ 
els, Fig. 2(b). In this case, and for a general orien¬ 
tation of iT, the Zeeman levels are the linear combina¬ 
tions of t and l spin states. As a result, transitions from 
each of the initial states on site i to both final states on 
site j are allowed. This fact distinguishes absorption of 
magnons from the conventional absorption of an ac elec¬ 
tric field 32-34 , and, as we will see below, gives rise to 
iT-dependence of the spin current. Another origin of in¬ 
dependence is the possibility of intrasite absorption of 
magnons at the boundary. We will see that the intrasite 
transitions dominate the absorption near the resonant 
condition fvjj = A z . Away from this condition, the inter¬ 
site transitions dominate. 


B. The model 


Consider a pair of localized states, i and j, Fig. 2(b). 
Assume for simplicity that the ferromagnet is an insu¬ 
lator, i.e. it is a barrier for electrons in N. Precession, 
m(t), of magnetization in ferromagnet can be modeled 
as a time-dependent correction oc m(t)& to the barrier 
potential. The pumping takes place since the wave func¬ 
tion, \Eq, can penetrate under the barrier. As a result, 
the Hamiltonian of site i has a correction 


SHi = J 


a x m x sin ujt + 6' y m y cos uot 


( 2 ) 


A. General considerations 

Figure 2 illustrates the difference between pumping 
into a metal, and into an insulator in an applied mag¬ 
netic field, H. While H is responsible for the magne¬ 
tization precession precession in the ferromagnet, it also 
causes a spin splitting, A 2 , of the spectrum in the metal¬ 
lic normal layer, Fig. 2a. This splitting, however, does 
not affect the absorption of magnons. The reason is that 
the absorption at a boundary does not require momen¬ 
tum conservation, i.e. the matrix element is constant, 
and thus there is no dependence of the spin current, I^ s \ 
on the dc field in the normal layer. 

The situation is different for an insulator, where the 


where J accounts for tunneling. Projections m x (t) and 
m y (t ) are proportional to the magnitude of the mi¬ 
crowave field and depend in a resonant way on the prox¬ 
imity of uj to uj FM r • Analytical expressions for these pro¬ 
jections can be found e.g. in Ref. [9]. 

The Hamiltonian SHi of Eq. (2) causes transitions of 
electrons between the sites i and j. Absorption of energy 
in course of these transitions is quite similar to the ab¬ 
sorption of the ac electric field by pairs of the localized 
states. However, the transitions caused by SHi are ac¬ 
companied by spin flips, both from j" to and from | to 
t- With regard to absorption of energy, one should add 
up the contributions of the both types of transitions, i.e. 

I {e) = h-> t + I^i (3) 
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However, the spin current results from the fact that these 
contributions are not equal to each other, so that 


I (s) = 4-t - (4) 

Thus, for calculation of the spin current into hopping 
insulator, one can use the standard “resonant” phonon¬ 
less absorption theory 32 and substitute the corresponding 
rates into Eq. (4). 


C. Resonant absorption at H = 0 

We first neglect the Zeeman splitting in the normal 
layer. In this case resonant transitions happen within 
pairs of localized states, Fig. lb. The correction SHi 
causes such transitions between the sites i and j because 
the corresponding wave functions | i) and | j) have a finite 
overlap integral, Uj 33 Due to this overlap, the eigenfunc¬ 
tions of the pair get modified as 

i»=i*>+ ^ 

|j> = ( 5 ) 

for Se = Ej — Si >0. The corresponding energies are 




T £ 


3 =r _ r = 
+ 2’ 


5s 2 + 4:t 2 j 


1 1/2 


(6) 


Since both modified eigenfunctions contain | i), the ma¬ 
trix element of 5Hi between them is finite, and the 
Golden-rule expression for the spin-flip part of the i j 
transition rate for ej > Si reads 


with lower energy (and there is a backflow of spins j. con¬ 
verted from t)- Since the occupation of the state at the 
interface is larger than of those at higher energy, there is 
a negative conversion rate because of transitions up 
the energy. This is exactly what Eq. (7) suggests. Fur¬ 
ther, since the occupation of the state at the interface 
is lower than of those at lower energy, there is a positive 
t conversion rate, or, again, negative one. Hence 
a simple permutation f(ii) f(sj) suffices to describe 
transitions to states with Sj < 

The product m x m y in Eq. (7) is specific for spin pump¬ 
ing, see Eq. (1). The expression for the net absorption 
rate contains ^(ra 2 + m 2 ) instead. Another difference 
from the conventional resonance absorption 31,32 is the 
structure of the matrix element in Eq. (7). This, how¬ 
ever, modifies the result of averaging over the sites, j, 
only by a numerical factor. A crucial observation in the 
averaging procedure 32 is that the relevant sites, j, are 
located within a narrow spherical layer with a radius 
which is found from the condition 2\tij(r UJ )\ = Huj. As¬ 
suming the exponential decay of the overlap integral with 
distance, \Uj(r)\ = to exp(—/a), we have 

i 2t ° rco 

r„=«ln-. (9) 

The result of averaging and summing over sites far away 
from the boundary reads 

/ (s) (w) = 27 T 2 m x m y J 2 ( guarl ) (10) 

where g is the density of states. The transition rate of 
Eq. (10) should be interpreted as the spin current gener¬ 
ated per a localized state coupled to the ferromagnet. 


D. Resonant absorption at finite H 


Ii% = -m x m y J 2 F(ii,ij,uj), (7) 

where the function F is defined as 


F {£ U £j,w) 


2 % HmtJM 

(ij — ii ) 2 — e i — hco ) 2 + (^) 2 

£[/&)-/&)] 

r 2 (r-M 2 + (f) 2 ' 


Here we have introduced the phonon broadening of the 
levels, t -1 . 

It is easy to see that the transition rate to states with 
Ej < Si is given by Eq. (7) with function F from Eq. (8), 
but with f{ii ) /(£j), and thus the rate has the same 

sign as Eq. (7). Physically, this can be seen from the fol¬ 
lowing argument: Consider the simple case of m x = m y . 
The Hamiltonian of Eq. (2) implies that for a given site 
at the interface, spins t are transferred to states of higher 
energy (and there is a backflow of spins t converted from 
| from those states), while spins 4- are pushed to states 


To generalize Eq. (7) to a finite magnetic field in the 
normal layer, one must take into account the modification 
of the spin eigenstates, as well as the Zeeman splitting in 
energies of the latter. The spin structure of the spin-split 
levels depends on the orientation of H as follows 

IXh+) = cos(^--- J|Xm+) +zsm(^--- 

(11) 

IXh-) = cos(^--- J|Xm-) +«sin^--- J|Xm+). 

( 12 ) 

Here the quantization axes for |xm±) and |xh±) spinors 
are chosen along the static part of the magnetization, and 
the external magnetic field, respectively, see Fig. 3(a). 
The states |xm±) at sites i and j are split by A z . 

All four transitions between states with |xm±) spin 
wave functions, Fig. 2(b), are allowed for a general ori¬ 
entation of the magnetic field. For spin-conserving tran¬ 
sitions (H-» + and- > —), the frequency dependence 
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of /remains cjr 2 , i.e. the same as in Eq. (10). Ori¬ 
entation of H enters into the prefactor: The product 
m x m y should be replaced with ^ sin 2 (Oh — 6 m )tu 2 f° r 
both transitions. 

While the spin-conserving transitions do affect the spin 
current density distribution in the sample, they are non¬ 
resonant, and it is the spin-flipping ones (+ —) that 

are responsible for the spin current generation at the in¬ 
terface. In other words, no spin current is possible in a 
stationary state without the latter processes. Therefore, 
in what follows we concentrate on the frequency and mag¬ 
netic field dependence of the corresponding rates. 

As far as + — and — + transitions are con¬ 
cerned, only the H-» — with absorption of a magnon, 

and-> + with emission of a magnon become important 

in the vicinity of the resonance hw = A z . The other two 
transitions are non-resonant, and therefore disregarded 
here. For the + <-> — transitions, the prefactor uj in the 
spin current remains intact, since it comes from the dif¬ 
ference in the populations of levels involved. However, 
despite the upper and lower Zeeman levels being sepa¬ 
rated in energy, the overlap of the spatial wave functions 
is determined by £*, ej in zero magnetic field. Thus, 
the + —>> — transitions take place between pairs with 
{sj — Si) ~ \1tlj — A z \. These pairs have the “shoulder” 


'I'huj — Az 


= a In 


2*q 

\hu - A z \' 


(13) 


Logarithmic divergence of Eq. (13), which is cut off at 
\fwj — A z \ ~ h/r, ensures the resonant character of spin¬ 
flipping transitions that we took into account. 

In addition to the replacement of r u by r u - a z in the 
spin current, the prefactor m x m y should be modified as 
m x m y —>> G(m x ,m y ), where the function G is defined as 

G(m x ,m y ) =\{m x + m y cos(0 H - 0m)) 2 , (14) 

so that the absorption, and thus the FMR damping, do 
not have the usual form oc m x m y . 

The most spectacular manifestation of the resonance 
hio = is that the intrasite transitions become pos¬ 

sible, as illustrated in Fig. 2(b). For these transitions 
the overlap of the spatial parts of the on-site wave func¬ 
tions is equal to 1, and the magnetic-field dependence of 
absorption is a pure Lorentzian. Orientation-dependent 
prefactor, which is the matrix element of SHi between 
the spinors \Xn+) and |% H _) is the same as in Eq. (14). 
Summarizing, we present the expression for spin current 
close to the resonance fuj = in the form 


J«(w) = 2G (m x ,m y )J 2 u^ 


J.Ji 2 


(a,-m 2 +(^ 


+ 7r V In 


\huj - A z 


, (15) 


where the first term comes from intrasite and the second 
term from intersite transitions. Directly at the resonance, 
the first term dominates. This is ensured by the condition 



FIG. 3. (Color online) (a) The geometry of FMR; For a fixed 
dimensionless frequency, Cj — cj F mr/4iyM S: the dimensionless 
magnitude, H = H/4tyM s , and orientation, Oh , of dc mag¬ 
netic field are related via Eq. (17). This dependencies are 
shown for the values of (D/y: (b) 0.5, (c) 1.2, and (d) 2. 

Red dots indicate the values of H , for which the condition 
7 H s = c o FM r is satisfied. 

(a) /(*)( u) [arb.units] (&) 




FIG. 4. (Color online) (a) the resonant condition yif = uj F mr 
is satisfied along the solid lines on the ) — Oh plane. The 

cutoff values of Oh are cos -1 ~ 55° and ty — cos 

(b) the behavior of the spin current calculated from Eq. (15) 
for ga 3 h/r = 4 ■ 10 -3 and tor/h = 15. 



ga 3 h/r 1. Since the combination 1 /go? is the minimal 
energy spacing between two sites in the insulator located 
within ^ a from each other, the above condition implies 
that this spacing is much bigger than the phonon broad¬ 
ening of individual levels, which is the definition of the 
Anderson insulator. As the deviation from the resonance 
increases, the behavior of I ^ (cj) is dominated by the sec¬ 
ond term. Neglecting the logarithm, the crossover takes 
place at \A Z — hcj\r/h > ( r/hga 3 ) 1. The behavior 

of spin current near the resonance is shown in Fig. 4(b), 
where the logarithm was cut off at \hw — A z \ = to/15. 


III. RESONANT ORIENTATIONS OF 
EXTERNAL FIELD 

Equation (15) is our main result. To make connec¬ 
tion to the experimental papers Refs. [23-26], below we 
calculate the magnitude and orientation of the dc field 
where the anomalous behavior of ISHE voltage takes 
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place. Such behavior takes place when two conditions 
are met: The Zeeman splitting of the localized states is 
equal to hoo, and uo = oo FMR . 

We specify the orientation of H and magnetization, 
M, using the notations common in the literature, see 
e.g. Refs. [9, 15, and 24], and Fig. 3. We will also intro¬ 
duce dimensionless variables H , M and Co, which stand 
for H, M and uo FMR in the units of 47 tM s , where M s is 
the saturation magnetization. Then the angle Om , corre¬ 
sponding the equilibrium orientation of M, is found from 
the condition that M is parallel to the effective magnetic 
field, with the demagnetizing term taken into account 9 


Pk 



FIG. 5. (Color online) “Spin-resistor” network. Polarizations 
Pi and Pj on the sites i, j determine the spin current be¬ 
tween these sites. The coefficients, Rij , are proportional to 
the electric hopping resistances. 


2 H sin {0 H - 0 M ) + sin 2 0 M = 0 , (16) 

while the expression for the resonant frequency, 4), 
reads 28 


<Co\ 2 _ 
K'y J 


H cos (Oh — 0 m) ~ cos 2 0m 

H cos (0 H - 0 M ) ~ cos 2 0 M • (17) 


From these two equations we exclude 0m and plot the di¬ 
mensionless field H versus Oh , for a given FMR frequency 
Co. Examples of these curves are shown in Fig. 3. Res¬ 
onant orientation is obtained by crossing a curve H(Qh) 
by the line Co = jH. Two intersections determine the ori¬ 
entations for which uo FMR is equal to the Zeeman splitting 
of the localized states. Upon changing uo FMR , we get two 
lines of resonances, Fig. 4(a). They occupy two domains: 
0 < Oh < cos -1 ^ and ( 7 r — cos -1 < Oh < tt. At 

the boundaries of the domains H goes to infinity. Then 
it follows from Eqs. (16) and (IT) that at these bound¬ 
aries sin (0m ~ Oh) = 0, and Oh satisfies the equation 

cos(20h) +cos 2 Oh = 0 , yielding Oh = cos -1 (^ 75 ) ~ 55°. 

In Refs. [23 and 24] on pumping into organics the mi¬ 
crowave frequency driving the resonance was 9.45 Ghz, 
while the values of 47 tM s were very different, namely, 
4 t tM s = 0.175 T in Ref. [23] and 4 t tM s = 0.805 T in 
Ref. [24]. Then from Fig. 4(a) we find that the resonant 
angle Oh should be close to 45° for Ref. [23] and 23° for 
Ref. [24]. 


IV. SPIN-RESISTOR NETWORK 

After the spin polarization is generated at the bound¬ 
ary, it should spread into the bulk of the insulator to 
avoid the backflow . 35 In a metal, where P is a contin¬ 
uous function of coordinates, this spreading is by spin 
diffusion accompanied by the Larmor precession. In a 
hopping insulator P takes discrete values, Pi, which are 
the polarizations on the sites, i. The Larmor precession 
is accounted for by the on-site Zeeman splitting, A z , of 
the levels, see Fig. 2. The frequencies of electron hops 


between two sites, i and j, depend strongly on their en¬ 
ergies, Si, Sj , and their spatial separation, r^. Then the 
issue of spreading of the spin polarization reduces to the 
question: what is the spin current • between the sites 
with polarizations Pi and Pj , provided that, on average, 
there is no charge current between these sites? 

If bias were applied between the two sites, then the 
average charge current, proportional to this bias, could 
be found by ascribing an effective resistance, , to the 
pair of sites 30,31 . It is easy to see that the same Re¬ 
determines the proportionality coefficient between l\ s _\- 
and Pi — Pj , namely 


rP) _ o ^ ^3 

± i^j ~ Z p dl ' 

^3 ds 


(18) 


In Eq. (18) we have assumed that the difference (^ — £j) 
is much smaller than the temperature, so that ^ is the 
same for both sites. Equation (18) follows from the fact 
that the on-site chemical potentials of the local majority 
and minority electrons are shifted by T 1-^1/§ 7 , respec¬ 
tively. The spinors that correspond to these local spin 
eigenstates are defined by the directions of Pi, Pj. Im¬ 
portantly, the fact that the chemical potential splitting is 
symmetric around the chemical potential of the unpolar¬ 
ized system ensures the absence of the charge current, i.e. 
the net current flow z —^ j is compensated by the net cur¬ 
rent flow j i. With different spin polarizations of the 
sites, the compensation of the charge flows leads to the 
imbalance of the spin flows, and thus to Eq. (18). Note 
that Eq. (18) remains valid in external magnetic field, 
which enters only via the magnitudes of polarizations. 
Overall, Eq. (18) suggests that polarization built up at 
the F-N boundary spreads along the current-carrying re¬ 
sistor network, as illustrated in Fig. 5. 


V. CONCLUDING REMARKS 

(i) Our result Eq. (15) applies when the phonon-induces 
broadening of the levels is smaller than uo. In the oppo¬ 
site case, uot 1, the mechanism of absorption is the 
Poliak-Geballe relaxation mechanism, Ref. [36]; no sharp 
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dependence of pumping near the resonance is expected 
in this regime. Unlike pumping into metals, the pumping 
rate Eq. (15) is not simply proportional to m x The 
prefactor in Eq. (14) depends on the relative orientation 
of m and the external magnetic field, H. 

(ii) Suppose that we are at resonance two = A z . The 
microwave field acts both by driving the FMR but also 
directly, by causing transitions between the Zeeman lev¬ 
els. If the amplitude of the field in frequency units (Rabi 
frequency) exceeds the inverse spin relaxation time, these 
transitions will be saturated in the bulk. Then the pump¬ 
ing becomes inefficient. 

(iii) In conventional theory of hopping transport the 
applied voltage drops not on all the resistors consti¬ 
tuting the network, but on the highest, critical, resis¬ 
tors representing the “hardest” hops. 31 The spin relax¬ 
ation rate will be dominated by hyperfine or spin-orbit 


environment 37,38 of this hop. 

(iv) We did not consider effects of electron-electron in¬ 
teraction, and did not describe in detail how finite resis¬ 
tance of the spin-current network, Section IV, affects the 
measured value of the spin current. In brief, Coulomb 
correlations enhance the absorption of magnons by in¬ 
creasing the number of singly occupied pairs 33 , while the 
measured spin current is given by Eq. (15) only in the 
limit of vanishingly small bulk resistance. These consid¬ 
erations are, however, completely standard, and do not 
change the qualitative picture of spin current generation 
by resonant magnon absorption in hopping insulators, 
developed in this paper. 
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